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Abstract: Total Internal Reflection Microscopy (TIRM) is a sensitive 
non-invasive technique to measure the interaction potentials between a 
colloidal particle and a wall with femtonewton resolution. The equilibrium 
distribution of the particle-wall separation distance z is sampled monitoring 
the intensity / scattered by the Brownian particle under evanescent illumi- 
nation. Central to the data analysis is the knowledge of the relation between 
/ and the corresponding z, which typically must be known a priori. This 
poses considerable constraints to the experimental conditions where TIRM 
can be applied (short penetration depth of the evanescent wave, transparent 
surfaces). Here, we introduce a method to experimentally determine I{z) 
by relying only on the distance-dependent particle-wall hydrodynamic 
interactions. We demonstrate that this method largely extends the range of 
conditions accessible with TIRM, and even allows measurements on highly 
reflecting gold surfaces where multiple reflections lead to a complex I{z). 
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1. Introduction 

Total Internal Reflection Microscopy (TIRM) |[T1 |2l is a fairly new technique to optically mea- 
sure the interactions between a single colloidal particle and a surface using evanescent light 
scattering. The distribution of the separation distances sampled by the particle's Brownian mo- 
tion is used to obtain the potential energy profile U{z) of the particle-surface interactions with 



suh-ksT resolution, where ksT is the thermal energy. Amongst various techniques available 
to probe the mechanical properties of microsystems, the strength of TIRM lies in its sensitiv- 
ity to very weak interactions. Atomic Force Microscopy (AFM) JSlIU requires a macroscopic 
cantilever as a probe and is typically limited to forces down to several piconewton 
the sensitivity of Photonic Force Microscopy (PFM) 15] |6l |71 can even reach a few femtonew- 
tons but this method is usually applied to bulk measurements far from any surface. 

TIRM, instead, can measure forces with femtonewton resolution acting on a particle near a sur- 
face. Over the last years, TIRM has been successfully appUed to study electrostatic lHIH, van 
der Waals ifTol . depletion ifTTlfTSlfTslfT?! . magnetic flSl . and, rather recently, critical Casimir 
lfT6l forces. 
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Fig. 1. Total Internal Reflection Microscopy (TIRM). (a) Schematic of a typical TIRM 
setup: a Brownian particle moves in the evanescent electromagnetic field generated by to- 
tal internal reflection of a laser beam; its scattering is collected by an objective lens; and 
the scattering intensity is recorded using a photomultiplier (PMT). (b) Typical experimental 
scattering intensity time-series (polystyrene particle in water, R = 1.45 flm). (c) Exponen- 
tial intensity-distance relation (/3 = llOnin). (d) Particle position distribution (acquisition 
time 1200i, sampling rate 500 /fz). (e) Experimental (dots) and theoretical (line) potential 
obtained from the position distribution using the Boltzmann factor. 

A schematic sketch of a typical TIRM setup is presented in Fig.lTJa). To track the Brownian 
trajectory of a spherical colloidal particle diffusing near a wall, an evanescent field is created 
at the substrate-liquid interface. The scattered light is collected with a microscope objective 
and recorded with a photomultiplier connected to a data acquisition system. Fig. [TJb) shows a 
typical example of an experimentally measured intensity time-series /, of a polystyrene particle 
with radius R— 1.45 lim in water 

Due to the evanescent illumination, the intensity of the light scattered by the particle is quite 
sensitive to the particle-wall distance. If the corresponding intensity-distance relation I{z) is 
known (and monotonic), the vertical component of the particles trajectory zt can be deduced 
from To obtain I{z), it is in principle required to solve a rather complex Mie scattering 
problem, i.e. the scattering of a micron-sized colloidal particle under evanescent illumination 



close to a surface ifTTlfTSl . where multiple reflections between the particle and the substrate and 
Mie resonances must be accounted for |[T9ll20l 12111221 . When such effects can be neglected, the 
scattering intensity is proportional to the evanescent field intensity and, since the latter decays 
exponentially, TIRM data are typically analyzed using a purely exponential I{z) — loe^'l ^ H] 

|2l[T7l[l8l|23l (Fig. [Tic)), where j3 = A /ATZ^n^sin^ B -n^ is the evanescent field penetration 
depth, A the incident light wavelength, the substrate refractive index, «„, the liquid medium 
refractive index, and 9 the incidence angle, which must be larger than the critical angle 9c =■ 
arcsin(«,„/n5). /q is the scattering intensity at the wall, which can be determined e.g. using a 
hydrodynamic method proposed |[24| . 

From the obtained z, the particle-wall interaction potential, i.e. is easily derived by 

applying the Boltzmann factor = —kBT Inp(z) to the calculated position distribution p{z) 
(Fig.[nd,e)). For an electrically charged dielectric particle suspended in a solvent, the interac- 
tion potential typically corresponds to t/(z) = Bexp(— K"z) + [|7r/?^(pp — p,„)^ — F^z. The first 
term is due to double layer forces with JC^' the Debye length and B a prefactor depending on 
the surface charge densities of the particle and the wall |[I]|2]|9l- The second term describes the 
effective gravitational contributions with Pp and p„, the particle and solvent density and g the 
gravitational acceleration constant; Fs takes into account additional optical forces, which may 
result from a vertically incident laser beam often employed as a two-dimensional optical trap 
to reduce the lateral motion of the particle 1251 . Depending on the experimental conditions, 
additional interactions, such as depletion or van der Waals forces, may arise. 

Despite the broad range of phenomena that have successfully been addressed with TIRM, 
most studies have been carried out with small penetration depths (at most )3 w lOOnm), and 
have therefore been limited to rather small particle-substrate distances z. In addition, no TIRM 
studies on highly reflecting walls, e.g. gold surfaces, have been reported, although such surfaces 
are interesting since they can support surface plasmons enhancing the evanescent field f26| and 
the optical near-field radiation forces 1271 l28l |29l . Furthermore, gold coatings can be easily 
functionalized l30l . which would allow to apply TIRM to e.g. biological systems. The reason 
for these limitations is the aforementioned problem to obtain a reliable l{z) relationship under 
these conditions. For example, it has been demonstrated that large penetration depths (e.g. 
above « 200 «m in Ref. l3T| ) increase the multiple optical reflections between the particle and 
the wall, which in turn leads to a non-exponential /(z) 1211 [3TI . Experiments combining TIRM 
and AFM found deviations from simple exponential behaviour very close to the wall even 
for shorter penetration depths 1321 . In principle, such effects can be included into elaborate 
scattering models, however, this requires precise knowledge of the system properties and, in 
particular, of the refractive indices of particle, wall, and liquid medium ISTll . Since the latter 
are prone to significant uncertainties (in particular for the colloidal particles), the application of 
TIRM under such conditions remains inaccurate. 

Here, we introduce a method to experimentally determine l{z) by making solely use of the 
experimentally acquired /, and of the distance-dependent hydrodynamic interactions between 
the particle and the wall. In particular, no knowledge about the shape of the potential IJ {£) is 
required. We demonstrate the capability of this method by experiments and simulations, and we 
also apply it to experimental conditions with long penetration depths (j3 = 720«m) and even 
with highly reflective gold surfaces. 

2. Theory 

2.1. Diffusion coefficient and skewness ofBrownian motion near a wall 

Colloidal particles immersed in a solvent undergo Brownian motion due to collisions with sol- 
vent molecules. This erratic motion leads to particle diffusion with the Stokes-Einstein diffusion 
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Fig. 2. Vertical diffusion coefficient £'x(z) near a wall (Eq. Q). 



coefficient D^e ~ hgT /6Kr]R, where Tj is the shear viscosity of the liquid. It is well known that 
this bulk diffusion coefficient decreases close to a wall due to hydrodynamic interactions. From 
the solution of the creeping flow equations for a spherical particle in motion near a wall assum- 
ing nonslip boundary conditions and negligible inertial effects, one obtains for the diffusion 
coefficient in the vertical direction ||33]| . 

= (1) 

l{z) 



2sinh((2«+l)a(z))+(2n+l)sinh(2a(z)) _ , 
4 sintf ( (n+0.5) a(z))- (2n+ 1 )2 sinh- (a(z)) 



and 



where l{z) = | sinh ( a (z) ) , ^^n-itiil^+i) 
a{z) ~ cosh^' + shown in Fig. |2l first increases with z approaching the 

corresponding bulk value at a distance of several particle radii away from the wall. 

Experimentally, the diffusion coefficient can be obtained from the mean square displacement 
(MSD) calculated from a particle trajectory. For the z-component this reads {{zt+ht ^ Zt)^) ~ 
2DsE^t, where (...) indicates average over time t. To account for a z-dependent diffusion coef- 
ficient close to a wall, one has to calculate the conditional MSD given that the particle is at time 
t at position z, i.e. ((zr+A/ —Zt)^ \zt = z) ~ 2Dx(z)Af where the equality is valid for Af 0; 
in such limit, this expression is only determined by the particle diffusion even if the particle is 
exposed to an external potential U{z)- From this follows that £>x(z) can be directly obtained 
from the particle's trajectory 

Eq. (|2]i was employed already by several groups ||34l[35l to validate Eq. 

The distribution of particle displacements h (z; Zo , Af ) around a given distance zo converges to 
a gaussian for Af ^ and therefore its skewness - i.e. the normalized third central moment - 
converges to zero. Accordingly, 

-L/f"^- V I... Uo. (3, 

where M(z, Af) = {zt+Kt —Zt\zt=^z)= argmax /z(z;z, Af), where argmax indicates the argument 
that maximize the given function. 



2.2. Mean square displacement and skewness of the scattering intensity 



In a TIRM experiment, /i(z;zo,Af) is translated into a corresponding intensity distribution 
/!(/;/(), Af) around intensity Iq = I{zq), whose shape strongly depends on I{z). In Fig. [3] we 
demonstrate how a particle displacement distribution /2(z;zo,Af), which is gaussian for small 
Af, translates into the corresponding scattered intensity distribution /i(/;/o,A/) for an arbitrary 
non-exponential /(z) dependence. In the linear regions of I{z), the corresponding li{I;I(),At) are 
also gaussian with the half-width determined by the slope of the I{z) curve. In the non-linear 
part of I{z), however, a non-gaussian intensity histogram with finite skewness is obtained. 




Fig. 3. Relation between position distributions and intensity distributions. Brownian diffu- 
sion of a particle around a point is symmetric, leading for small At to a gaussian distribu- 
tion h(z',zo,At) (bottom). According to I{z) this leads to the scattered intensity histograms 
h{I;lQ,At) (left): in the linear region of I(z), h{I;I(,,At) is also gaussian with the width 
depending on /' (Eq. (|4j); in the nonlinear region of I{z) (central histogram), h{I;Io,At) 
deviates from a gaussian and has a finite skewness depending on /" (Eq. ([5j). 

In the following we calculate the MSD and the skewness of !i{I;Io,At) for an arbitrary I{z), 
which we assume to be a continuous function with well defined first and second derivates /' 
and /". In the vicinity of zq, I{z) can be therefore expanded in a Taylor series I{z) = I{zo) + 
I' (zo )(z-Zo) + \l" (zo )(z-Zo)^foTz^Zo, where /' = | and /" = || ■ The MSD of It (/; /q , Af ) 
is 

^ ' ' " (4) 



MSD(/) ^ lim ((/,^^, !/,=/)= I'^z) ■ 2Dx(z), 

At^O At 

where ( (1,+^ -I,Y\I,=l)= I'^ ( {zt+M -Zt)^ \I,= l) for Af ^ and Eq. dU has been used. 
The skewness of h{I;Io,At) is 



1 I,+Ar-Ir-M{I,At) 



aF-^q Af ^ \ \^ ^MSD(/ 
with M(/,Af) = argmax h{I;I,At] 



\It=I)=llMi-V2DKzj, (5) 



2|/'(z)| 

and ((/,+A,-/,-M(/,Af))3 |/,=/) 
/ 

I'mz,+A,-Zr-M{z,At))^\z,=z) + \l"-I"{{z,+Ar-Zt-M{z,At)f\zt^z) for Af 



where the first term is null because of Eq. (O, and the second term is calculated using the 
properties of the momenta of a gaussian distribution {{-f) = 3 ((...)^) ■ 

In Fig.|4]we applied Eqs. (|4| and Q to the intensity time-series corresponding to a particle 
trajectory simulated using a Langevin difference equation assuming various I{z) (Fig.|4ja,c,e)). 
The solid lines in Fig.|4l;b,d,f) show the theoretical MSD(/) (black) and S(/) (red) and the dots 
the ones obtained from the simulations. When I{z) is linear (Fig.Ua)), MSD(/) is proportional 
to Eq. ([T]i and S(/) vanishes (Fig.HJb)). When I{z) is exponential (Fig.HJc)) or a sinusoidally 
modulated exponential (Fig. SJe), MSD(/) is not proportional to Eq. ([T]i and large values of 
the skewness occur as shown in (Figs.|4td,f)). Small deviations between the theoretical curves 
and the numerical data can be observed for intensities where the particle drift becomes large in 
comparison to the time-step (At = 2ms); in our case this corresponds to a slope of the potential 
of about 1 pN/fim, which is close to the upper force limit of typical TIRM measurements. If 
necessary such deviations can be reduced employing shorter time-steps. 




Fig. 4. Various intensity-distance relations and their effect on MSD(/) (black) and skew- 
ness S(/) (red). Both theoretical values (solid lines) and the results from the analysis of 
numerically simulated data (dots) using Eq. Q and Eq. l|5j are presented (R = 1.45 /im, 
Pp = l.053g/cm^ , samples 10^, frequency lOOHz). (a)-(b) Linear /(z). (c)-(d) Exponen- 
tial /(z). (e)-(f) Exponential I{z) modulated by a sinusoidal function. 



2.3. Obtaining I{z) from It 

The correct I{z) satisfies the conditions 

r MSD(/(z)) = I'Hz)-2D^{z) 

where MSD(/) and S(/) are calculated from an experimental If. Thus, the problem of deter- 
mining I{z) can be regarded as a functional optimization problem, where Eqs. ^ have to be 
fulfiUed. 



3. Analysis workflow 

Here, we present a concrete analysis workflow to obtain z, from the experimental /, by finding 
the I{z) that satisfies Eqs. (|6]l. To do so, we will construct a series of approximations /''^(z) 
indexed by / converging to I{z). 

(1) As first guess, take /'^'(z) =/oexp(— z/j3)+^.v, where j3, Iq and bs are parameters chosen 
to optimize Eqs. Often some initial estimates are available from the experimental condi- 
tions: j3 can be taken as the evanescent field penetration depth, Iq as the scattering intensity at 
the wall, and bs as the background scattering in the absence of the Brownian particle. While j3 
is typically well known, Iq and b^ are prone to large experimental systematic errors and uncer- 
tainties. 

(2) Take /(i)(z) =/''''(z)[l -G(/('''(2),M*'',C7(i),A(i))], where G(x,^,(7,A) =Aexp(-(x- 
/x)^/a^) is a gaussian, and the parameters /i'^', a''', and A^'^ optimize Eqs. (|6]l. Gaussian 
functions were chosen because they have smooth derivatives and quickly tend to zero at infinite. 
Notice that the choice of a Gaussian is unessential for the working of the algorithm. 

(3) Reiterate step (2), substituting with and with until Eqs. Q are satisfied 
within the required precision. 

(4) Invert (z), i.e. numerically construct z^'^^^ (/). 

(5) Takez, =z('■+')(^^)■ 
4. Experimental case studies 

4.1. Validation of the technique 

We test our approach on experimental data (polystyrene particle with R=l.45 fXm near a glass- 
water interface kept in place by a vertically incident laser beam ||25]| ) for which the exponential 
/(z) is justified (j3 = 120«ot, X = 658 nm) [3 1 .1 . As illustrated in Fig. HJa), there is indeed 
agreement between the measured (dots) and theoretical potential (solid line). In the inset, the 
measured diffusion coefficient (black dots) agrees with Eq. ^ (black solid line) and the skew- 
ness (red dots) is negligible (small deviations in the region where the potential is steepest are 
due to the finite time-step). The criteria for /(z) in Eqs. ^ are already fulfilled after Iq and bs 
have been optimized in the first step of the analysis workflow in the previous section. As shown 
in Fig.|5lb), the experimental MSD(/) (black dots) and skewness S(/) (red dots) fit Eqs. (|4]i and 
(|5]i (solid lines). 

4.2. TIRM with large penetration depth 

We now apply our technique under conditions where an exponential /(z) is not valid, i.e. for 
large penetration depth as mentioned above. Fig. |6] shows the potential obtained for the same 
particle as in Fig.|5]but for a penetration depth (j3 = 720«m). Note, that compared to|5]the 
potential extends over a much larger distance range since the particle's motion can be tracked 
from hundreds of nanometers to microns. The green data points show the faulty interaction 
potential that is obtained when assuming an exponential /(z). Since the only difference is in 
the illumination, the same potential as in Fig. |5] should be retrieved (solid line in Fig. |6ta)). 
However, applying an exponential /(z) (green line in Fig. |6tb)) wiggles appear in the potential 
(green dots in Fig. |6la)). Their origin is due to multiple reflections between the particle and 
the wall as discussed in detail in ||2TI . The correct /(z) (black line in Fig. ISb)) is obtained 
with the algorithm proposed in the previous section: after 9 iterations the conditions in Eqs. (|6]l 
appear reasonably satisfied, as shown in Fig.|6jc). With this /(z), we reconstructed the potential 
represented by the black dots in Fig. |6ja), in good agreement with the one in Fig. |5] It should 
be noticed that, even though the deviations of the correct /(z) from an exponential function are 
quite small, this is enough to significantly alter the measurement of the potential. This again 



0.1 0.2 0.3 0.4 




0.25 0.5 0.75 1 



I{V\ 

Fig. 5. TIRM with exponential intensity-distance relation, (a) The experimental (dots) and 
theoretical (solid line) potential. Inset: the diffusion coefficient on the position data (black 
dots) fits well Eq. (T} (black solid line), while the absolute value of the Brownian motion 
skewness S(z) is small (red dots), (b) Experimental MSD(/) (black dots) and skewness 
S(/) (red dots) for a scattering intensity time-series (polystyrene particle, R = 1.45/xm, 
Pp = 1.053 g / cm^ , Hp = 1.59 suspended in water n„, = 1.33 with 300 /iM NaCl background 
electrolyte, ff^' = 17«m, near a glass surface ;?j = 1.52, acquisition time 1800^, sampling 
frequency 500Hz) calculated using Eq. ^ and Eq. (O. Given the short penetration depth 
(jS = 120;7m), the theoretical MSD(/) (black solid line) and S(/) (red solid line) for an 
exponential I{z) fit the experimental ones and the conditions in Eqs. ^ are fulfilled. 

demonstrates the importance of obtaining the correct I{z) for the analysis of TIRM experiments. 
4.3. TIRM in front of a reflective surface 

To demonstrate that our method is capable of correcting even more severe optical distortions, we 
performed measurements in front of a reflecting surface (20 «m gold-layer, reflectivity w 60%, 
j3 = 244 «m). The experimental conditions are similar to the previous experiments. Only the salt 
concentration was lowered to avoid sticking of the particle to the gold surface due to van der 
Waals forces, leading to a larger electrostatic particle-surface repulsion, and the optical trap was 
not used. Using an exponential /(z) (green line in Fig.|2tb)), we obtain the potential represented 
by the green dots in Fig.|7ta), which clearly features unphysical artifacts, e.g. spurious potential 
minima. After 27 iterations of the data analysis algorithm, the black I{z) in Fig.|7];b) is obtained, 
which reasonably satisfies the criteria in Eqs. ^ (Fig.lTjc)). The reconstructed potential (black 
dots in Fig. [Tja)) fits well to theoretical predictions (solid line); in particular the unphysical 
minima disappear. 

5. Conclusions & Outlook 

TIRM is a technique which allows one to measure the interaction potentials between a colloidal 
particle and a wall with femtonewton resolution. So far, its applicability has been limited by the 




I[V] 



Fig. 6. TIRM with large penetration depth, (a) The experimental potential (black dots) ob- 
tained using the fitted intensity-distance relation and the theoretical one (black solid line). 
The green dots represent the faulty potential obtained using the exponential /(z). (b) The 
fitted /(z) (black line) and the exponential one (green line) corresponding to the penetration 
depth j6 = llQnm. (c) Experimental intensity MSD(/) (black dots) and skewness S(/) (red 
dots) for a scattering intensity time-series (same particle and acquisition parameters as in 
Fig-HJ calculated using Eq. ^ and Eq. ([5}. Due to the large penetration depth, the I{z) 
diverges from an exponential; the theoretical MSD(/) (black solid line) and S(/) (red solid 
line) correspond to the fitted l{z)- 



need for an a priori knowledge of the intensity-distance relation. I{z) °= exp(— z/j3) can safely 
be assumed only for short penetration depths of the evanescent field and transparent surfaces. 
This, however, poses considerable constraints to the experimental conditions and the range of 
forces where TIRM can be applied. Here, we have proposed a technique to determine I{z) 
that relies only on the hydrodynamic particle-surface interaction (Eq. ^ and, differently from 
existing data evaluation schemes, makes no assumption on the functional form of I{z) or on 
the wall-particle potential. This technique will particularly be beneficial for the extension of 
TIRM to new domains. Here, we have demonstrated TIRM with a very large penetration depth, 
which allows one to bridge the gap between surface measurements and bulk measurements, 
and TIRM in front of a reflecting (gold-coated) surface, which allows plasmonic and biological 
applications. 

This new technique only assumes the knowledge of the particle radius, which is usually 
known within an high accuracy and can also be measured in situ ll24ll . and the monotonicity of 
/(z). Were I{z) not monotonous, as it might happen for a metallic particle in front of a reflective 
surface, the technique can be adapted to use the information from two non-monotonous signals, 
e.g. the scattering from two evanescent fields with different wavelength liTTl . We notice that the 
technique encounters its natural limits when Eq. ([T]) does not correctly describe the particle-wall 
hydrodynamic interactions. This may happen in situations when the stick boundary conditions 
do not apply or when the hydrodynamic interactions are otherwise altered, e.g. in a viscoelastic 
fluid. 
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Fig. 7. TIRM in front of a reflective surface, (a) The experimental potential (black dots) 
obtained using the fitted intensity-distance relation and theoretical one (solid black line). 
The green dots represent the faulty potential obtained using the exponential I{z) with 
j8 = 244nm. (b) The fitted I{z) (black line) and the exponential one (green line) corre- 
sponding to the evanescent field penetration depth j8 = 244 nm. (c) Experimental intensity 
MSD(/) (black dots) and skewness S(/) (red dots) for a scattering intensity time-series 
(same particle and acquisition parameters as in Fig. [S] except for background electrolyte 
50 fxM, K^^ = 42«m) calculated using Eq. ^ and Eq. (|5]l. Due to the presence of a 
20;!m-thick gold layer on the surface, the I(z) deviates from an exponential; the theoretical 
MSD(/) (black solid line) and S(/) (red solid line) correspond to the fitted I{z). 



Since the conditions in Eqs. (|6|l are fulfilled only by the correct I{z), they permit a self- 
consistency check on the data analysis. Even when an exponential I{z) is justified, errors that 
arise from the estimation of some parameters (e.g. the zero-intensity Iq and the background 
intensity bs) can be easily avoided by checking the consistency of the analyzed data with the 
aforementioned criteria. In principle, the analysis of TIRM data can be completely automatized, 
possibly providing the missing link for a widespread application of TIRM to fields, such as 
biology, where automated analysis techniques are highly appreciated. 

The proposed technique can also be useful to determine the intensity-distance relation in 
all those situations where it is possible to rely on the knowledge of the system hydrodynam- 
ics, while the scattering is not accurately known. Often explicit formulas are available for the 
hydrodynamic interaction of an over-damped Brownian particle in a simple geometry, while 
complex numerical calculations are needed to determine its scattering. As a limiting case, this 
technique might also prove useful for the PFM technique working in bulk, where the diffusiv- 
ity is constant. Indeed, under certain experimental conditions - e.g. using back-scattered light 
instead of the more usual forward-scattered light ll36l - the intensity-distance relation can be 
non-trivial and it can be necessary to determine it experimentally. 



